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ABSTRACT 

CR7 is the brightest z = 6.6 Lya emitter (LAE) known to date, and spectroscopic 
follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III 
stars. We examine this interpretation using cosmological hydrodynamical simulations. 
Several simulated galaxies show the same “Pop III wave” pattern observed in CR7. 
However, to reproduce the extreme CR7 Lya/HeII1640 line luminosities (La/Heii) a 
top-heavy IMF and a massive ( ^ IO^Mq) Pop III burst with age 2 Myr are re¬ 
quired. Assuming that the observed properties of Lya and Hell emission are typical 
for Pop HI, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs 
at z = 6.6 with La > lO'^^'^ergs”^ should also host Pop HI stars producing an ob¬ 
servable LHeii ~ 10^^'^erg s“^. As an alternate explanation, we explore the possibility 
that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). 
Our model predicts La, LHeii, and X-ray luminosities that are in agreement with the 
observations. In any case, the observed properties of CR7 indicate that this galaxy 
is most likely powered by sources formed from pristine gas. We propose that further 
X-ray observations can distinguish between the two above scenarios. 
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1 INTRODUCTION 

The end of the Dark Ages is marked by the appearance of the 
first stars. Such - Pop III - stars had to form out of a pris¬ 
tine composition (H-(-He) gas with virtually no heavy ele¬ 
ments. Lacking these cooling agents, the collapse had to rely 
on the inefficient radiative losses provided by H 2 molecules. 
Mini-halos, i.e. nonlinear dark matter structures with mass 
Mh ~ 10®“^Mq collapsing at high redshift (z ~ 30), are 
now thought to be the preferred sites of first star formation 
episodes (Yoshida et al. 2006; Turk et al. 2009; Salvador! & 
Ferrara 2009; Greif et al. 2012; Visbal et al. 2015). Although 
the Initial Mass Function (IMF) of Pop III stars is largely 
uncertain, physical arguments suggest that they could have 
been more massive than present-day (Pop II) stars. Further¬ 
more, the metals produced by Pop III stars polluted the sur¬ 
rounding gas (Bromm et al. 2002; Wise et al. 2012; Xu et al. 
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2013), inducing a transition to the Pop II star formation 
mode (“chemical feedback”, Schneider et al. 2002, 2006). 
Metal enrichment is far from being homogeneous, and pock¬ 
ets of pristine gas sustaining Pop III star formation can in 
principle persist down to z ~ 3 — 4 (Tornatore et al. 2007; 
Trent! et al. 2009; Maio et al. 2010; Salvador! et al. 2014; 
Pallottini et al. 2014a; Ma et al. 2015), yielding Pop III star 
formation rate (SFR) densities of ~ lO“"^M 0 yr“^Mpc“®, i.e. 

1% of the Pop II SFR density at those redshifts. 

The search effort for Pop III stars at moderate and high 
redshifts has become increasingly intense in the last few 
years (e.g. Kashikawa et al. 2012; Heap et al. 2015). Observa- 
tionally, a galaxy hosting a recent (t* 2 Myr) Pop III star 

formation episode should show strong Lya and He ii fines 
and no metal lines (e.g. Schaerer 2002; Raiter et al. 2010; 
Kehrig et al. 2015). Until now, no indisputable evidence for 
Pop III stars in distant galaxies has been obtained, and ob¬ 
servations have only yielded upper bounds on Pop HI SFR 
(e.g. Cai et al. 2011; Cassata et al. 2013; Zabl et al. 2015). 
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This situation might dramatically change following the re¬ 
cent observations of CR7 by Sobral et al. (2015, S15 here¬ 
after). 

CR7 is the brightest Lya emitter (LAE) at 2 > 6 , and 
it is found in the COSMOS field (Matthee et al. 2015). Spec¬ 
troscopic follow-up by S15 suggests that CR7 might host a 
PopIII-like stellar population. This is based on the aston¬ 
ishingly bright Lya and He ii lines {La — 10^®'®®ergs“^, 
Lueii — 10‘*^'^®ergs“^) and no detection of metal lines. S15 
shows that CR7 can be described by a composite of a PopIII- 
like and a more normal stellar population, which would have 
to be physically separated, and that would be consistent 
with e.g. Tornatore et al. (2007). HST imaging shows that 
CR7 is indeed composed of different components: 3 sepa¬ 
rate sub-systems (A, B, C) with projected separations of 
^ 5kpc. FllOW(YJ) and F160W(H) band photometry in¬ 
dicates that clump A might be composed of young (blue) 
stars, while the stellar populations of B-l-C are old and 
relatively red. The observed Lya and He ii lines are nar¬ 
row (FWHM < 200kms"^ and FWHM < 130kms“S re¬ 
spectively), disfavoring the presence of an AGN or Wolf- 
Rayet (WR) stars, which are expected to produce much 
broader (FWHM .^5 10®kms“^) lines (e.g. De Breuck et al. 
2000; Brinchmann et al. 2008; Erb et al. 2010). S15 con¬ 
cluded that CR7 likely contains a composite stellar popu¬ 
lation, with clump A being powered by a recent Pop III- 
like burst (t* 2Myr), and clumps B-l-C containing an old 
(t* ~ 350Myr) burst of Pop II stars with M* ~ lO^'^M©, 
largely dominating the stellar mass of the entire system. 

Based on cosmological simulations that follows the si¬ 
multaneous evolution of Pop II and Pop III stars (Pallottini 
et al. 2014a, P14 hereafter), we examine the interpretation 
of CR7 as a Pop HI host system and explore its implications. 
We also propose an alternate explanation, briefly discussed 
in S15, where CR7 is powered by accretion onto a Direct 
Collapse Black Hole and suggest further tests. 


2 SIMULATION OVERVIEW 

We use the ACDM cosmological^ hydrodynamical simula¬ 
tions presented in P14 (see that paper for a comprehen¬ 
sive description), obtained with a customized version of the 
Adaptive Mesh Refinement code RAMSES (Teyssier 2002) 
to evolve a (10/i“^ Mpc)® volume from 2 = 99 to 2 = 4, 
with a dark matter mass resolution of ~ 5 x lO®/i“^M 0 , 
and an adaptive baryon spatial resolution ranging from 
~ 20 / 1 “^ kpc to ~ lh~^ kpc. Star formation is included 
via sub-grid prescriptions based on a local density thresh¬ 
old. If the star forming cell gas has metallicity below (above) 
the critical metallicity, Zcrit = 10“"^Z©, we label the newly 
formed stars as Pop HI (Pop II). Supernova feedback ac¬ 
counts for metal-dependent stellar yields and return frac¬ 
tions appropriate for the relevant stellar population^. The 

^ We assume a ACDM cosmology with total matter, vacuum 
and baryonic densities in units of the critical density Ha = 
0.727, ^dm = 0.228, f!;, = 0.045, Hubble constant Hq = 
100 h km s“^ Mpc“^ with h = 0.704, spectral index n = 0.967, 
(Tg = 0.811 (Larson et al. 2011). 

^ While in P14 we explore different types of IMF for Pop III, here 
we show results assuming a Pop Il-like Salpeter IMF. It is to note 



Figure 1. Metallicity (Z) map centered on the simulated galaxy 
MB45. Pop II locations are marked with filled stars, for recent 
(blue, f* 2 Myr) and old (red) formation events, respectively. 
Locations of the recent (old) Pop III burst are shown by the open 
blue (red) pentagons. Overdensity (log A = 1, 2, and 3) contours 
are plotted with black lines. The scale is indicated in arcsecond 
and kpc. 

simulated galaxy sample reproduces the observed cosmic 
SFR (Bouwens et al. 2012; Zheng et al. 2012) and stellar 
mass density (Gonzalez et al. 2011) evolution in the red- 
shift range 4 sC 2 10, and - as shown in Pallottini et al. 

(2015) - P14 reproduce the observed the luminosity func¬ 
tion at 2 = 6 . Additionally, the derived Pop HI cosmic SFR 
density is consistent with current observational upper limits 
(e.g. Nagao et al. 2008; Cai et al. 2011; Cassata et al. 2013). 
To allow a direct comparison with CR7 we will concentrate 
on the analysis of the 2 ~ 6 simulation output. 

2.1 Pop Ill-hosting galaxies 

As noted in S15, the interpretation of CR7 fits in the “Pop 
III wave” scenario suggested by Tornatore et al. (2007). As 
an example of Pop III wave in action, we show the case 
of “MB45”, a simulated P14 galaxy with total stellar mass 
M* = lO^'^M©. In Fig. 1, we plot the metallicity {Z) and 
overdensity (A) map around MB45. The star formation his¬ 
tory in MB45 starts with a Pop III event. These stars ex¬ 
plode as supernovae enriching with metals the central re¬ 
gions of MB45. As a result, star formation there continues in 
the Pop II mode, while in the less dense external regions, not 
yet reached by the metal-bearing shocks. Pop HI stars can 
still form. The process repeats until the unpolluted regions 
have densities exceedingly low to sustain star formation. 

The total (i.e. old-|-young stars) Pop HI mass in MB45 
is M 3 ~ 1 O® ®M 0 ; about 20 % of this stellar mass formed 
in a recent burst (age ^2 Myr). The total stellar mass 
(M* ~ 10®M©) of MB45 is dominated by Pop II stars pro¬ 
duced at a rate SFR 2 — O. 5 M 0 yr“^. Thus, while MB45 

that our simulations suggests that the Pop III SFR seems almost 
independent from the IMF (see in particular Fig. 14 in P14). 
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formation activity proceeds in the Pop III wave mode and 
resembles that of CR7, the physical properties of MB45 and 
CR7 are different, because of the 2 orders of magnitude dif¬ 
ference in total stellar mass. A direct comparison between 
the PopIII-PopII separation in CR7 (projected ~ 5 kpc) and 
MB45 (10 kpc, projected 5 kpc), although fairly consis¬ 
tent, might not be very meaningful due to the different mass 
of the two systems. This is because the separation depends 
on the mass-dependent metallicity profile in galaxy groups 
(see Fig.l in Pallottini et al. 2014b). However, we note that 
our current simulated volume is simply to small to be able to 
directly recover sources such as CR7, with volume densities 
of ~ 10"®Mpc"^ 

The simulated volume contains many galaxies hosting 
Pop III star formation, whose number density (rigai) as a 
function of their total stellar mass (M*) is shown in Fig. 2, 
along with their Pop III mass (M 3 ). Pop III stars are found 
preferentially in low-mass systems M* IO^'^Mq, that typ¬ 
ically form these stars in a series of a few M 3 ~ 1 O®M 0 
bursts^, before Pop III formation is quenched by chemical 
feedback. 

Note that all larger galaxies, M* IO^Mq, contain 
some Pop III component inherited from progenitor halos. 
We can then regard M 3 calculated by summing old and 
young Pop III stars as a solid upper bound to the total 
Pop III mass produced during the galaxy lifetime. The size 
of the simulation volume, dictated by the need of resolv¬ 
ing the very first star-forming halos, is too small to fairly 
sample the mass function of galaxies with M* L 10®Mq and 
more massive. To make predictions in this high-mass range, 
we slightly extrapolate the simulated trends of Ugai and M 3 
(lines in Fig. 2). We caution that such extrapolation might 
imply uncertainties. 


2.2 Lya and He ii emission 

We label as pure ( composite) galaxies whose emission is pro¬ 
duced by Pop II only (Pop II-|-Pop III) stars. For a pure 
galaxy, the Lya line luminosity is (e.g. Dayal et al. 2008): 

= 2.8 X 10^^ SFR 2 ergs"^, (la) 

where Aa is an attenuation factor that accounts for both 
internal (interstellar medium) absorption and intergalactic 
medium transmissivity; Pop II SFR is expressed in MQyr“^ 
and to relate this to the M* we assume SFR 2 /M* = sSFR ~ 
2.5Gyr“^, consistent with our simulations and observations 
(Daddi et al. 2007; McLure et al. 2011; Gonzalez et al. 2011, 
see Sec. 3). Pure galaxies have no He ii emission. For a com¬ 
posite galaxy, the Lya emission is the sum of the contribu¬ 
tions from Pop II and Pop III components: 

+ AJ^Ma , (lb) 

where I 3 is the Pop HI Lya line luminosity per unit stellar 
mass'^. Analogously, the He ii emission is given by 

Lneii = , (Ic) 

where is the Pop III He ii luminosity per unit stellar 
mass. Both I 3 and depend on the IMF and burst age, 

^ We refer to App. A of P14 for possible resolution effects. 

We are implicitly assuming that Pop III stars form in a burst, 
an assumption justified by the analysis presented in Sec. 2.1. 



log(/M./M.) 


Figure 2. Upper panel: Number density (rigai) of simulated 
galaxies as a function of their total stellar mass (AL* = M 2 + M 3 ) 
for: (i) all galaxies (pink bar), (ii) galaxies with old+young Pop 
III stars (red bars), (iii) galaxies with a young Pop III component 
(blue bars). All quantities are averaged on log(AAL/M 0 ) ~ 0.3 
bins. Black dot-dashed and dotted lines correspond to the best-fit 
analytical extrapolation of rtgai- Lower panel: Mass of Pop III 
{M 3 ) in composite galaxies considering young (blue circle) and 
old-l-young (red squares) stars. In both panels the vertical dashed 
line indicates the value of M* for CR7 as inferred by SED fitting 
(see Sec. 2.2). 


f*. We adopt the Pop III models by Schaerer (2002) and 
Raiter et al. (2010), and we use a Salpeter IMF (power-law 
slope a = —2.35), with variable lower {miow) and upper 
(mup) limits. As long as m„p L 1 O^M 0 , the results are very 
weakly dependent on the upper limit, which we therefore fix 
to rriup = 1 O®M 0 , leaving miow as the only free-parameter. 

As noted by S15, to reproduce CR7 La and Lueii with 
Pop His, a mass of M 3 ~ 1 O^“®M 0 newly-born {tir A 2 — 5 
Myr) stars is required, depending on the IMF. Such a large 
amount of young Pop HI stars is contained in none of P14 
galaxies and it is not predicted by the adopted analytical ex¬ 
trapolation (see Fig. 2). Thus none of the simulated compos¬ 
ite galaxies would reproduce CR7 line emission. However, it 
is possible that CR7 might have experienced a more vigorous 
Pop HI star formation burst as a result of a very rare event - 
e.g. a recent major merger - not frequent enough to be cap¬ 
tured in our limited box volume. As an estimate we adopt 
the value of M 3 resulting from the sum of all (old-l-young) 
Pop HI stars formed in our galaxies. 

Under this hypothesis, we can hx miow, by using the 
zero age main sequence (ZAMS) tracks (t* = 0). By SED 
fitting, S15 shows that CR7 Pop II stellar mass (completely 
contained in clumps B-l-G) is likely M* ~ M 2 ~ 1 O^°M 0 . 
From the lower panel of Fig. 2, we find that this mass 
corresponds to a Pop HI mass of M 3 ~ 1 O^'®M 0 . From 
Pop HI SED hts of region A S15 estimate M 3 ~ 1 O^M 0 . 
As these stars must be located in CR7 clump A, whose 
He II luminosity is Lueii = lO'^^'^ergs”^, eq. Ic requires that 
^Heii _ lo^^ ^ergs”^ Mg^. In turn this entails a top-heavy 
IMF with miow = 6.7 M 0 . 

Having fixed the IMF, we can readily derive the pre¬ 
dicted Pop HI contribution to the Lya emission; this turns 
out to be I 3 = 10^®'^ergs“^ Mg^. GR7 has an observed 
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Figure 3. Upper panel: Number of LAEs (Na) as a function 
of their Lyo luminosity {La) in the COSMOS/UDS/SA22 fields 
(Matthee et al. 2015, green pentagons). The histogram shows the 
predicted number of composite (i.e. Pop III hosts) LAEs in the 
same survey; the bar colors indicate the expected Lnell emission, 
also shown by the numbers in the bars. Lower panel: Fraction 
of composite galaxies eq. 2 a) for which the correction for 

the finite volume bias (eq. 2b) of Matthee et al. (2015) survey has 
been included (red solid) or neglected (black dashed). 


La = 10^®'®ergs“^, with no contribution from clumps B+C 
(SFRg = 0). This comparison allows us to determine, using 
eq. lb, the Lya line attenuation factor, Aa = 10“°'®^. 

Roughly 66% of the line luminosity is therefore damped, 
a figure consistent with other derivations (e.g. Dayal et al. 
2008), and with the analysis of S15. The above procedure 
provides a basis to model Lya and He ii emission for both 
pure and composite galaxies, assuming that the properties 
of the Pop III component are similar to those derived from 
CR7. 


served volume®. Therefore, it is important to account for 
the statistical (Poisson) fluctuations of the galaxy number 
counts as follows: 

1463(1 ± , (2b) 

where k = (composite, pure). The distribution of ncomp(iH*) 
is shown in the upper panel of Fig. 2 as the “young” Pop 
III curve®; while Ppme accounts for the remaining galaxies. 
The effect of the finite volume effects on can be ap¬ 

preciated from the lower panel of Fig. 3. Assuming a higher 
sSFR = 5Gyr“^ (e.g. Stark et al. 2013) yields the modifica¬ 
tions shown by the dash-dotted line. 

In the upper panel of Fig. 3 we plot the LAE number 
{No) as a function of La- Matthee et al. (2015) observations 
(green pentagons) are shown along with our predictions for 
the composite LAE and expected LhbII emission. CR7 is 
the most luminous LAE observed, and it is in the brightest 
luminosity bin {La = 10"^"^^®'®ergs“^); by assumption CR7 
is a composite galaxy. If so, we then predict that out of the 
46 (30) LAEs with La = 10'‘®'^±° ®ergs"^ (> 10^® ®ergs“S 
cumulative), 13 (14) must also be composite galaxies^, with 
observable LhbII — 10"'^'®ergs“^ ( 4 10'‘^'^ergs“'^). Follow¬ 
up spectroscopy of those luminous Lyman-a emitters at 2 ; = 
6.6 will allow to test this prediction. We recall that this test 
assumes that all Pop III give raise to the same Lya and 
He II emission as inferred from CR7, that requires that all 
the Pop HI stellar mass was formed in a single burst with 
age 4 2 Myr. 

Particularly in the regime where /comp < 1 (see lower 
panel in Fig. 3), a sample of LAEs is needed to test 
our model predictions. Eor example for “Himiko”, the sec¬ 
ond most luminous® confirmed LAE at 2 = 6.6 with 
La — 10"^® "^ergs“^ (Ouchi et al. 2009), for which recent 
VLT/X-Shooter observations have provided a 3cr limit of 
Lneii 4 lO^^'^ergs”^ (Zabl et al. 2015), our model predicts 
LhbII — lO'^^'^ergs”^ i.e. a four times higher He ii luminos¬ 
ity. However, this is predicted only for /comp — 20 — 30% of 
galaxies at this La- 


3 PREDICTIONS FOR BRIGHT LAES 

Starting from the assumption that CR7 is a “typical” com¬ 
posite galaxy, and using mio-w = 6.7 M© and Aa = 10“®'®^, 
we can now predict how many LAEs among those ob¬ 
served by Matthee et al. (2015) are composite galaxies, 
i.e. contain Pop HI stars. The number of LAEs in the 
COSMOS/UDS/SA22 fields with luminosity La is Na = 
^a{La) Voba, where ^a is the observed Lya luminosity func¬ 
tion and 1463 = 4.26 x 10® Mpc® is the observed volume. 
Among these, a fraction 

f.omp ^ Naon.p/{N,amp + Apme) ( 2 a) 

contain Pop HI stars, where Acomp (Npure) is the number of 
composite (pure) galaxies in 146s at a given La- 

Eq.s la and lb, show that a given La can be produced 
by a composite galaxy with a lower M* with respect to 
a pure (PopH) galaxy. Eor instance. La — 10‘^® ®ergs“^ 
requires M* ~ 1 O^® ®M 0 for a pure galaxy, but only 
M* 4 1 O®'®M 0 for a composite one. Such large objects are 
very rare at the redshift of CR7 (2 = 6 . 6 ) and in the ob- 


4 ALTERNATIVE INTERPRETATION 

Given the extreme conditions required to explain the ob¬ 
served properties of GR7 in terms of Pop HI stars and a 
set of assumptions, it is worth exploring alternative inter¬ 
pretations. The most appealing one involves Direct Gollapse 
Black Holes (DGBH), which is briefly discussed in S15. High- 
2 pristine, atomic halos {Mh 4 1 O®M 0 ) primarily cool via 
Lya line emission. In the presence of an intense Lyman- 
Werner (LW, E-y = 11.2— 13.6 eV) irradiation, Hg molecule 
photo-dissociation enforces an isothermal collapse (Shang 

® As a reference, a M* ~ 1O^® ®M 0 is hosted in a dark mat¬ 
ter halo of mass ~ 1 O^^'®M 0 , whose abundance is nj, ~ 
10“®Mpc“® at 2 ~ 6 (e.g. Sheth & Tormen 1999). 

® As noted in Sec. 2.1, all galaxies with M* > 10^ have old Pop 
III stars, thus considering the “old+young” track for ricomp would 
yield an unrealistically high composite number. 

^ The predicted number would become 7 (7), by assuming sSFR 
from Stark et al. (2013). 

® As re-computed in S15 using Y band to estimate the continuum, 
in order to match the calculation for CR7. 
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et al. 2010; Latif et al. 2013; Agarwal et al. 2013; Yue et al. 
2014), finally leading to the formation of a DCBH of initial 
mass M, ~ 1O‘* ®“® ®M 0 (Begelman et al. 2006; Volonteri 
et al. 2008; Ferrara et al. 2014), eventually growing np to 
1 O®“^M 0 by accretion of the halo leftover gas. 

In CR7, clump A appears to be pristine, and it 
is irradiated by a LW flux from B+C® of ~ 5 x 
10 “^®ergs“^ cm“^Hz“^ sr“^, well in excess of the required 
threshold for DCBH formation (Shang et al. 2010; Latif et al. 
2013; Sugimura et al. 2014; Regan et al. 2014). Thus CR7 
might be a perfect host for a DCBH. 

We investigate the time-evolving spectrum of an accret¬ 
ing DCBH of initial mass M, = 1 O®M 0 by coupling a ID 
radiation-hydrodynamic code (Pacucci & Ferrara 2015) to 
the spectral synthesis code CLOUDY (Ferland et al. 2013), as 
detailed in Pacucci et al. (2015). The DCBH intrinsic spec¬ 
trum is taken from Yue et al. (2013). The DCBH is at the 
center of a halo of total gas mass Mg ~ 1 O^M 0 , distributed 
with a core plus a r~^ density profile spanning up to 10 
pc. The accretion is followed until complete depletion of the 
halo gas, i.e. for ~ 120 Myr. During this period the total 
absorbing column density of the gas varies from an initial 
value of ~ 3.5 x 10^"^cm“^ to a final value ^ 10^^cm“^, i.e. 
from mildly Compton-thick to strongly Compton-thin. Note 
that while Lya attenuation by the interstellar medium is in¬ 
cluded, we do not account for the likely sub-dominant IGM 
analogous effect. 

Fig. 4 shows the time evolution of the Lya, He ii and 
X-ray (0.5-2 keV) luminosities. Both Lya and He ii are con¬ 
sistent with the observed CR7 values during an evolutionary 
phase lasting ~ 17 Myr (14% of the system lifetime), longer 
than the shorter period (t* L 2 Myr) of our assumption for 
a massive Pop HI burst. 

The equivalent width of the He ii line in the CR7 com¬ 
patibility region ranges from 75 to 85 A. The column density 
during the CR7-compatible period is ~ 1.7 x lO^'^cm”^, i.e. 
mildly Compton-thick. The associated X-ray luminosity is 
A lO^^ergs”'^, fully consistent with the current upper limit 
for CR7 (L lO'^'^ergs"^, Elvis et al. 2009). Deeper X-ray 
observations of CR7 might then confirm the presence of the 
DCBH. However, this limit is already obtained with 180 ks 
of integration time on Chandra, meaning that a stringent 
test might only be possible with the next generation of X- 
ray telescopes. 


5 CONCLUSIONS 

CR7 is the brightest 2 = 6.6 LAE in the COSMOS field 
(Matthee et al. 2015). Spectroscopic follow-up (Sobral et al. 
2015) suggests that CR7 might host Pop HI stars, along 
with Pop II and thus be explained by a “Pop HI wave” 
scenario. We have further investigated such interpretation 
using cosmological simulations following the formation of 
Pop II and Pop HI stars in early galaxies. 

We find simulated galaxies (like MB45 in Fig. 1) host¬ 
ing both Pop III and Pop II stars at 2 = 6.0. Such “compos¬ 
ite” galaxies have morphologies similar to that of CR7 and 

® The LW is estimated by accounting for the stellar properties of 
clump B-l-C (in particular see Fig. 8 in S15), and by assuming a 
5 kpc distance between B-l-C and A. 



Figure 4. Upper panel: Time evolution of the Ly« (blue 
solid line). He II (red dashed line) line and X-ray (violet dot 
dashed line) luminosities calculated for the accretion process onto 
a DCBH of initial mass 1O®M0. The green shaded region indi¬ 
cates the period of time during which our simulations are compat¬ 
ible with CR7 observations. The current upper limit for X-ray is 
% 10'^'^ ergs“^, (Elvis et al. 2009, horizontal violet line). Lower 
panel: Time evolution of the Hell/Lya lines ratio. The black 
horizontal dashed line indicates the observed values for CR7. 

consistent with the “Pop HI wave scenario”. However, to re¬ 
produce the extreme CR7 Lya/HeHI640 line luminosities, 
a top-heavy IMF combined with a massive (M 3 L 1 O^M 0 ) 
Pop HI burst of young stars L 2 — 5 Myr) is required. 
Our simulations do not predict such large burst, i.e. M 3 ~ 
1 O®M 0 , but our volume is also smaller than that used to dis¬ 
cover CR7. Nonetheless, assuming that CR7 is typical of all 
metal-free components in our simulations, we predict that 
in the combined COSMOS, UDS and SA22 fields, out of the 
30 LAEs with La > lO^^'^ergs”^, 14 should also host Pop 
HI stars producing an observable I/Heii ~ 10'*^'^ergs“^. 

Given the extreme requirements set by the Pop HI in¬ 
terpretation, we explored the possibility that CR7 is instead 
powered by accretion onto a Direct Collapse Black Hole 
(DCBH) of initial mass 1 O®M 0 . The predicted La and I/ueii 
match CR7 observations during a time interval of ~ 17 Myr 
(~ 14% of the system lifetime). The predicted CR7 lumi¬ 
nosity at 0.5-2 keV, L lO^^ergs”^, is significantly below 
the current upper limit, i.e. 10 "^^ergs“^. 

We conclude that the DCBH interpretation of CR7 is 
very appealing, and competitive with the explanation in¬ 
volving a massive Pop HI burst. For both explanations, the 
dominant ionizing source of this galaxy should have formed 
from pristine gas. Deep X-ray observations and other follow¬ 
up observations should allow to shed more light on this very 
peculiar source. 


ACKNOWLEDGMENTS 

SS acknowledges support from the Netherlands Orga¬ 
nization for Scientihc research (NWO), VENI grant 
639.041.233. RS acknowledges support from the European 
Research Council under the European Union (FP/2007- 


MNRAS 000, 1-6 (2015) 









6 Pallottini et al. 


2013)/ERC Grant Agreement n. 306476. D. Sobral ac¬ 
knowledges (i) financial support from the NWO through 
a Veni fellowship and (ii) funding from FCT through 
a FCT Investigator Starting Grant and Start-up Grant 
(IF/01154/2012/CP0189/CT0010) and from FCT grant 
PEst-OE/FIS/UI2751/2014. 


References 

Agarwal B., Davis A. J., Khochfar S., Natarajan P., Dunlop J. S., 
2013, MNRAS, 432, 3438 

Begelman M. C., Volonteri M., Rees M. J., 2006, MNRAS, 370, 
289 

Bouwens R. J., et al., 2012, ApJ, 754, 83 

Brinchmann J., Pettini M., Chariot S., 2008, MNRAS, 385, 769 
Bromm V., Coppi P. S., Larson R. B., 2002, ApJ, 564, 23 
Cai Z., et al., 2011, ApJL, 736, L28 
Cassata P., et al., 2013, A&A, 556, A68 
Daddi E., et al., 2007, ApJ, 670, 156 

Dayal P., Ferrara A., Gallerani S., 2008, MNRAS, 389, 1683 
De Breuck C., Rottgering H., Miley G., van Breugel W., Best P., 
2000, A&A, 362, 519 
Elvis M., et al., 2009, ApJS, 184, 158 

Erb D. K., Pettini M., Shapley A. E., Steidel C. C., Law D. R., 
Reddy N. A., 2010, ApJ, 719, 1168 
Ferland G. J., et al., 2013, Revista Mexicana de Astronomia y 
Astrofisica, 49, 137 

Ferrara A., Salvadori S., Yue B., Schleicher D. R. G., 2014, 
preprint, (arXiv:1406.6685) 

Gonzalez V., Labbe L, Bouwens R. J., Illingworth G., Franx M., 
Kriek M., 2011, ApJL, 735, L34 

Greif T. H., Bromm V., Clark P. C., Glover S. C. O., Smith R. J., 
Klessen R. S., Yoshida N., Springel V., 2012, MNRAS, 424, 
399 

Heap S., Bouret J.-C., Hubeny L, 2015, preprint, 
(arXiv:1504.02742) 

Kashikawa N., et al., 2012, ApJ, 761, 85 

Kehrig C., Vilchez J. M., Perez-Montero E., Iglesias-Paramo J., 
Brinchmann J., Kunth D., Durret F., Bayo F. M., 2015, ApJL, 
801, L28 

Larson D., et al., 2011, ApJS, 192, 16 

Latif M. A., Schleicher D. R. G., Schmidt W., Niemeyer J. C., 
2013, MNRAS, 436, 2989 

Ma Q., Maio U., Giardi B., Salvaterra R., 2015, MNRAS, 449, 
3006 

Maio U., Ciardi B., Dolag K., Tornatore L., Khochfar S., 2010, 
MNRAS, 407, 1003 

Matthee J., Sobral D., Santos S., Rottgering H., Darvish B., 
Mobasher B., 2015, MNRAS, 451, 4919 
McLure R. J., et al., 2011, MNRAS, 418, 2074 
Nagao T., et al., 2008, ApJ, 680, 100 
Ouchi M., et al., 2009, ApJ, 696, 1164 
Pacucci F., Ferrara A., 2015, MNRAS, 448, 104 
Pacucci F., Ferrara A., Volonteri M., Dubus G., 2015, preprint, 
(arXiv:1506.05299) 

Pallottini A., Ferrara A., Gallerani S., Salvadori S., D’Odorico V., 
2014a, MNRAS, 440, 2498 

Pallottini A., Gallerani S., Ferrara A., 2014b, MNRAS, 444, L105 
Pallottini A., Gallerani S., Ferrara A., Yue B., Vallini L., Maiolino 
R., Feruglio C., 2015, preprint, (arXiv: 1506.05803) 

Raiter A., Schaerer D., Fosbury R. A. E., 2010, A&A, 523, A64 
Regan J. A., Johansson P. H., Wise J. H., 2014, ApJ, 795, 137 
Salvadori S., Ferrara A., 2009, MNRAS, 395, L6 
Salvadori S., Tolstoy E., Ferrara A., Zaroubi S., 2014, MNRAS, 
437, L26 

Schaerer D., 2002, A&A, 382, 28 


Schneider R., Ferrara A., Natarajan P., Omukai K., 2002, ApJ, 
571, 30 

Schneider R., Omukai K., Inoue A. K., Ferrara A., 2006, MNRAS, 
369, 1437 

Shang C., Bryan G. L., Haiman Z., 2010, MNRAS, 402, 1249 
Sheth R. K., Tormen G., 1999, MNRAS, 308, 119 
Sobral D., Matthee J., Darvish B., Schaerer D., Mobasher B., 
Rottgering H., Santos S., Hemmati S., 2015, ApJ, 808, 139 
Stark D. P., Schenker M. A., Ellis R., Robertson B., McLure R., 
Dunlop J., 2013, ApJ, 763, 129 

Sugimura K., Omukai K., Inoue A. K., 2014, MNRAS, 445, 544 
Teyssier R., 2002, A&A, 385, 337 

Tornatore L., Ferrara A., Schneider R., 2007, MNRAS, 382, 945 
Trenti M., Stiavelli M., Michael Shull J., 2009, ApJ, 700, 1672 
Turk M. J., Abel T., O’Shea B., 2009, Science, 325, 601 
Visbal E., Haiman Z., Bryan G. L., 2015, preprint, 

(arXiv:1505.06359) 

Volonteri M., Lodato G., Natarajan P., 2008, MNRAS, 383, 1079 
Wise J. H., Turk M. J., Norman M. L., Abel T., 2012, ApJ, 745, 
50 

Xu H., Wise J. H., Norman M. L., 2013, ApJ, 773, 83 
Yoshida N., Omukai K., Hernquist L., Abel T., 2006, ApJ, 652, 6 
Yue B., Ferrara A., Salvaterra R., Xu Y., Chen X., 2013, MNRAS, 
433, 1556 

Yue B., Ferrara A., Salvaterra R., Xu Y., Chen X., 2014, MNRAS, 
440, 1263 

Zabl J., Nprgaard-Nielsen H. U., Fynbo J. P. U., Laursen P., 
Ouchi M., Kjaergaard P., 2015, MNRAS, 451, 2050 
Zheng W., et al., 2012, Nature, 489, 406 

This paper has been typeset from a Te^/IMT)X file prepared by 
the author. 


MNRAS 000, 1-6 (2015) 


